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We consider a magnetic s=l/2 impurity in the antiferromagnetic spin chain as a function of two 
coupling parameters: the symmetric coupling of the impurity to two sites in the chain Ji and the 
coupling between the two sites J2. By using field theory arguments and numerical calculations we 
can identify all possible fixed points and classify the renormalization flow between them, which leads 
to a non-trivial phase diagram. Depending on the detailed choice of the two (frustrating) coupling 
strengths, the stable phases correspond either to a decoupled spin with Curie law behavior or to a 
non-Fermi liquid fixed point with a logarithmically diverging impurity susceptibility as in the two 
channel Kondo efi'ect. Our results resolve a controversy about the renormalization fiow. 
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Impurities in low dimensional magnetic and elec- 
tronic systems tiave recently received a lot of atten- 
tion in the context of mesoscopic systems, quasi one- 
dimensional compounds and high temperature supercon- 
ductivity. Much progress has been made in particular 
for one-dimensional systems, where quantum impurities 
are known to renormalize as the temperature is lowered, 
so that the impurity can be described by an effective 
boundary condition in the low temperature limit |^ . The 
temperature dependence of the impurity contributions to 
the specific heat and the susceptibility can also be pre- 
dicted. Very often the antiferromagnetic spin-1/2 chain 
with impurities is taken as the prototype example for this 
renormalization behavior P-[l^, which captures the es- 
sential renormalization behavior of many other impurity 
systems. 

A number of different impurity configurations have 
been considered where typically one coupling parameter 
is varied. The chain with one altered link is known to 
renormalize to an effective open boundary conditions |^] 
in analogy with the renormalization of the conductivity 
in a quantum wire JTzt . However, two neighboring al- 
tered links renormalize to a periodic boundary condition 
1 2| in analogy with resonant tunneling in quantum wires 
1 18[ . This renormalization behavior is also equivalent to 
the two-channel Kondo effect [0^. Equally interesting 
is the coupling of an external impurity spin to one or two 
sites in the chain. An impurity spin s coupled to one 
site in the chain gets screened with a resulting decoupled 
singlet of spin s — 1/2 and an open chain with one site 
removed [^,^-^. An impurity spin coupled to two sites 
in the chain has first been considered in Ref. ||T^ where 
an equivalence with the two channel Kondo effect was 
found on the basis of the field theory description ||f^,|fl| , 
which would result in an "overscreened" impurity spin 
with a logarithmically diverging impurity susceptibility. 
However, numerical simulation studies later questioned 



this result and predicted a decoupled impurity spin as 
the stable boundary condition instead with a Curie law 
susceptibility ||^,^ , which is an unresolved controversy. 

We now consider a two parameter impurity model. 
This allows us to consider the different impurity config- 
urations above in a richer phase diagram with a greater 
number of possible fixed points. Using field theory tech- 
niques and advanced numerical results we are able to map 
out the phase separation line between an overscreened 
and a decoupled impurity spin exactly, which not only 
resolves the controversy, but also shows a more complex 
renormalization behavior. 

The Hamiltonian we consider here describes a Heisen- 
berg chain where two neighboring sites are coupled to an 
impurity spin-1/2 with strength Ji. The two sites are 
also coupled to each other with strength J2 as shown in 
Fig.0 

H = J Si ■ Si+i + JiSq ■ (^Sn + Sij + J2SN ■ Si. 
1=1 

(1) 
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FIG. f. Impurity model with two parameters Ji and J2. 

In the low temperature limit this system is known to 
be well described by by a level 1 Wess-Zumino-Witten 
(WZW) model with a marginal irrelevant operator [ p^ . 
A spin operator at position x in the chain can be ex- 
pressed in terms of the current operators J and the WZW 
field g 
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S{x) Jl + Jr + (-l)"" const. tYag. (2) 

General values of Ji and J2 break the conformal invari- 
ance of the effective field theory and introduce perturbing 
operators in the field theory Hamiltonian, which in turn 
can be classified by their renormalization behavior. As 
the temperature is lowered the system always must ap- 
proach a fixed point that restores conformal invariance 
[01 . For the spin-1/2 chain in Eq. (|l|) the possible fixed 
points are always given in terms of simple boundary con- 
ditions on the spin chain (open or periodic). In addi- 
tion there may be completely decoupled impurity spins 
present. The possible fixed points for this system are 
shown in Fig. |^. 

• Ji = J; J2 = 0. A periodic chain with + 1 sites 
and no impurity spin. The impurity spin has been 
absorbed. We denote this fixed point by Pn+i- 

• Ji = 0, J2 — 0. An open chain with N sites and 
a decoupled s=l/2 impurity spin. We denote this 
fixed point by On ® \- 

• Ji = 0, J2 — J- A periodic chain with N sites 
and a decoupled s=l/2 impurity spin, denoted by 

Pn®!- 

• Ji = 0, J2 = 00. An open chain with N — 2 sites 
and a decoupled impurity s=l/2. The two end sites 
at 1 and N have locked into a singlet state. We 
denote this fixed point by On-2 ® f • 




FIG. 2. The four possible fixed points. 

The renormalization behavior in the vicinity of each of 
those fixed points is determined by the operator content 
of the impurity Hamiltonian according to the perturba- 
tion of Ji and J2 around the corresponding fixed point 
values. Typically the operators with the lowest scaling 
dimensions (leading operators) can be determined by ei- 
ther a straight-forward symmetry analysis or by direct 
application of Eq. (||). Local operators become relevant 
when their scaling dimensions is less than one d < 1. 

In the vicinity of the periodic fixed point Pn+i the 
leading irrelevant operator is given by dx^T^g because site 
parity symmetry does not allow any more relevant oper- 
ators [p|. We therefore immediately conclude that this 



fixed point is stable in all directions in the Ji — J2 phase 
diagram, i.e. small perturbations on both Ji and J2 are 
irrelevant as indicated in the phase diagram Fig |[ 




FIG. 3. The full phase diagram of the Ji — J2 impurity 
model. The possible fixed points are indicated by the thick 
solid dots. The points on the phase separation line have been 
estimated numerically by analyzing {Si ■ Sn) as in Fig. ^ for 
132 difi'erent values for Ji and J2. 



In the vicinity of the the fixed point On ® ^ the fields 
in the perturbing Hamiltonian are given by boundary 
operators. This stems from the analytic continuation 
of the left moving fields in terms of right moving fields 
in order to restore the conformal invariance at an open 
boundary. The leading operator (Jl(0) -f- Jl{N)) ■ Si^ip 
of scaling dimension d = 1 is created by the coupling 
Ji to the impurity spin from the open ends. This op- 
erator is marginally relevant for antiferromagnetic sign 
and marginally irrelevant for ferromagnetic sign analo- 
gous to the Kondo effect ||^. The coupling between the 
end spins J2 can only produce the irrelevant operators 
II J^(0), J^(iV) and Ji(0) • Jl{,N) with scaling dimen- 
sion d — 2. Hence the coupling J2 is always irrelevant in 
the vicinity of On ® k- We therefore conclude that the 
fixed point On ® 5 is stable for Ji < 0, but for Ji > 
the effective coupling to the impurity spin increases un- 
der renormalization and flows towards the stable fixed 
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point Pn+1, corresponding to a "healing" of the chain 
M as indicated in the phase diagram Fig I This effect 
is in fact completely analogous to the renormalization 
flow and overscreening of the spin excitations in the two- 
channel Kondo problem resulting in a logarithmi- 
cally diverging impurity susceptibility. 

Near the fixed point On-2 (8> ^ the impurity spin is 
separated by a locked singlet in the limit J2 — > oo which is 
effectively decoupled from the rest of the chain (unless we 
introduce additional coupling constants). Therefore the 
fixed point appears to be stable in all directions 

in our Ji — J2 phase diagram in Fig ^ 

A more interesting scenario can be found near the fixed 
point Pn ® § since a small coupling Ji of either sign is 
frustrating. Here the most relevant operator trg corre- 
sponds to a slight modification of one link in the chain 
J2 . If the impurity spin is absent this operator is respon- 
sible for the breaking of the chain , but the frustrating 
coupling Ji must now also be considered. In particu- 
lar, a small coupling to the impurity spin Ji produces 
the operator {Jl + Jr) ■ Simp which is marginally rele- 
vant/irrelevant for antiferromagnetic/ferromagnetic sign 
respectively. The irrelevant operator dxtrag ■ S'lmp of di- 
mension d = 3/2 is also created by Ji. In summary we 
may express the impurity Hamiltonian at the fixed point 
Pn ® 5 as 

^^imp = 7iti-g + 72(/l + Jr) ■ ^imp + l3dj;tTag ■ S^mp, 



where to lowest order 



71 cx ( J2 - J), 72 oc Jl, 73 cx Jl. 



(3) 



(4) 



The renormalization group equations can then be de- 
termined by using the operator product expansion. We 
find that 
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(5) 
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where the dot indicates the derivative in respect to the 
logarithm of the cutoff. The renormalization flow is dom- 
inated by the relevant coupling constant 71, so that it is 
crucial to determine if 71 is initially negative or positive. 
Interestingly, the irrelevant coupling constant 73 plays an 
important role in defining the corresponding phase sepa- 
ration line along 71 — 37273, which together with Eq. (Q) 
gives a parabolic shape in the Ji — J2 phase diagram 
J2 — J oc J^ as indicated in Fig. ^. Indeed if the irrelevant 
coupling constant 73 is neglected, one might erroneously 
come to the conclusion that the relevant backscatter- 
ing 71 is exactly zero along the line J2 = J resulting 
in just a marginal renormalization fiow as postulated in 
Refs. [pOUll) which is not the case. 



If 71 > the effective coupling J2 increases quickly so 
that the two end spins lock into a singlet independent 
of the value of Ji. Therefore, the stable fixed point is 
On-2 ® 5 above the phase separation line. If, however, 
7i < the effective coupling J2 decreases towards zero 
and the stable fixed point now depends on the marginal 
coupling Jl. For Ji < the stable fixed point is On <8i ^ 
while for Ji > the "healing" process towards the fixed 
point Pn+1 takes place again. The leading irrelevant op- 
erators close to Pn+1 will again produce a logarithmically 
divergent impurity susceptibility , but the renormaliza- 
tion fiow is not completely analogous to the two channel 
Kondo effect since the operator content near the unstable 
fixed point is quite different. In particular the relevant 
operator trg is completely absent in the ordinary Kondo 
problem. The impurity susceptibility in the phases char- 
acterized by the fixed points On ® 5 and On-2 ® 5 is a 
Curie law behavior as T ^ from the decoupled impu- 
rity spin degrees of freedom. 

We have now determined the basic shape of the phase 
diagram, but to quantitatively map out the phase sepa- 
ration line it is necessary to resort to numerical methods. 
We therefore choose the transfer matrix renormalization 
group method for impurities, which can determine local 
correlation functions directly in the thermodynamic limit 
Q. One obvious order parameter is the correlation be- 
tween the two end spins which takes well-defined values 
at each fixed point in the zero temperature limit 
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Sn) = J 
Sn)=0 
Sn) = I -ln2 
Sn) = -J, 



41n2+|C(3) 



(6) 



where C is the Riemann Zeta function. The values at 
the periodic chain fixed points correspond to the first 
and second nearest neighbor correlations in the ordinary 
Heisenberg chain pO| ], while the open chain fixed point 
values correspond to two uncorrelated spins and a spin 
singlet, respectively. Numerical simulations can never 
reach zero temperature in the thermodynamic limit, but 
it is well possible to determine if the correlation between 
the spins renormalizes towards larger or smaller values 
as the temperature is lowered and thereby determine the 
stable fixed point. 

In particular we postulated that above the phase- 
transition line in Fig. ^ the two end spins effectively lock 
into a singlet so that the correlation should decrease with 
decreasing temperature. Below the phase transition line, 
however, the two spins will only be weakly correlated at 
On 5 or ferromagnetically correlated at Pn+i so that 
the correlation should increase with decreasing tempera- 
ture. 
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FIG. 4. The correlation {Si ■ Sn) as a function of T for 
Ji = 0.8J and J2 = 0, 0.4, 0.8, 1, 1.2, 1.5, 1.8, 2J from 
above. By analyzing the slope or absolute value compared 
to I — In 2 the phase transition can be estimated to be at 
J2 « 1.44 J in this case. 

This behavior is shown in Fig. |^ for Ji = 0.8 J and 
different values of J2. As T ^ the correlation clearly 
decreases for larger values of J2 ^ 1.5 J but increases for 
smaller values J2 ^ 1.2 J so that the phase transition 
occurs somewhere between those two values. We can 
identify the critical value of J2 more closely by consider- 
ing the slope and absolute value of {Si ■ Sn) as T — * 0. 
The critical value can be defined as the point where the 
slope is zero or the absolute value is 1/4 — In 2, which 
gives approximately the same result of J2 « 1.44 J. We 
have made a similar analysis for many different points 
in the Ji — J2 parameter space and thereby mapped out 
the exact numerical location of the phase separation line 
as shown in Fig. As Ji, J2 — * 00 the phase transition 
occurs at Ji — J2 which can be determined by analyzing 
the ground state of the three coupled spins Sn, Sq, Si. 

We now wish to relate our results to previous numer- 
ical studies of this impurity model by Zhang et al |^|-^ , 
where three points in the phase diagram were analyzed 
at (Ji = 0.4 J, J2 = J), (Ji = 0.4 J, J2 = 0), and 
(Ji = 1.5 J, J2 — J), which led to a controversy with 
the predictions in Refs. ||lo| , pT[| . In particular, the numer- 
ical studies found no logarithmically diverging impurity 
susceptibility at (Ji = 0.4 J, J2 = J). As mentioned 
above the analytic studies in Refs. failed to con- 

sider the leading relevant operator trg at the unstable 
fixed point Pn ® 5, but the logarithmic impurity sus- 
ceptibility is created by the leading irrelevant operator 
dxtrg near the stable fixed point Pn+i Q, so that we also 
postulate a logarithmic divergence below some crossover 
temperature Tk- However, looking at the phase diagram 
in Fig. H we see that the point (Ji = 0.4J, J2 = J) 
lies very close to the phase transition line, so that the 



crossover temperature Tr- must be extremely small and 
cannot be reached by any numerical method today, which 
explains the controversy. On the other hand the points 
(Ji = 1.5, J2 = J) and (Ji = 0.4, J2 = 0) are fur- 
ther from the phase separation and the crossover tem- 
perature is larger. Therefore, the logarithmic diver- 
gence can indeed be observed ^^j^ ■ In Ref. ||] the end 
spin correlation {Si ■ Sn) was also determined, which 
is consistent with our results, but with a different in- 
terpretation. Because the overall value of {Si ■ Sn) is 
negative it was concluded in that work that the point 
(Ji = 0.4J, J2 — J) produces a Curie law susceptibil- 
ity from frustration, i.e. an effectively decoupled impu- 
rity spin. However, the temperature dependence clearly 
shows a renormalization towards an increasing correla- 
tion as T ^ 0, so that Pn+i is the correct stable fixed 
point with an absorbed impurity spin and a logarithmi- 
cally diverging susceptibility below an extremely small 
crossover temperature Tk ■ Interestingly, there is an inte- 
grable spin-1 chain model with a spin-1/2 impurity which 
can also reproduce this two-channel Kondo renormaliza- 
tion behavior near the fixed point Pn+i |pl|| . However, 
the interesting boundary fixed points that were recently 
discovered in this integrable model have no analogous 
expression in our parameter space. 

In conclusion we have analyzed a complex impurity 
model in the spin-1/2 chain which shows a rich phase 
diagram. We were able to identify three different sta- 
ble phases that are attracted to well-defined fixed points 
under renormalization. The three phases are separated 
by the unstable fixed point Pn ® ^ which shows a high 
sensitivity to the detailed choice of couplings due to frus- 
tration. Using field theory arguments and numerical 
techniques we were able to map out the phase diagram 
exactly and also resolve a controversy in the literature 
|f7|,p pC|JTl[ | . Interestingly, an irrelevant operator plays a 
crucial rule in determining the shape of the phase sepa- 
ration line. It would be interesting to see if the current 
findings can be generalized to other systems where a frus- 
trated state may separate more stable phases (e.g. in a 
ladder system with a zigzag geometry) . In general, a frus- 
trated state should always be unstable since it has been 
postulated that renormalization typically occurs towards 
a state with a lower ground state degeneracy |2^. The 
system we considered here is a typical example of this 
behavior. 
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